Investigating the exclusive protoproduction of dileptons at high energies 
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I. INTRODUCTION 

Recently, the physics of large impact parameter inter- 
actions at the LHC and Tevatron has raised great inter- 
est (see Refs. fl] for a review). These electromagnetic 
interactions in pp and AA collisions extend the physics 
program of photon induced processes beyond the energies 
currently reached at DESY-HERA. At very high collision 
energies, the electromagnetic field surrounding a nucleus 
contains photons energetic enough to produce new par- 
ticles in ultra-peripheral collisions. This can happen in 
a purely electromagnetic process through a two-photon 
interactions or in an interaction between a photon from 
one of the nuclei and the other "target" nucleus. The 
high photon energies and fluxes lead to large cross sec- 
tions for several photon-induced reactions. The coherent 
contribution from the Z protons in the nucleus, further- 
more, enhances the number of equivalent photons by a 
factor Z^. 

These ultraperipheral collisions (UPCs) are a good 
place to constraint the photonuclear cross sections as 
the dominant processes in UPCs are photon-nucleon (nu- 
cleus) interactions. The current LHC detector configu- 
rations can explore small- a; hard phenomena with nuclei 
and nucleons at photon-nucleon center-of-mass energies 
above 1 TeV, extending the x range of HERA by a factor 
of ten. In particular, it will be possible to probe diffrac- 
tive and inclusive parton densities in nuclei using several 
processes. For instance, the interaction of small dipoles 
with protons and nuclei can be investigated in elastic and 
quasi-elastic J/tjj and T production as well as in high t 

production accompanied by a rapidity gap. Several of 
these phenomena provide clean signatures of the onset of 
the new high gluon density QCD regime [2]. The LHC is 
in the kinematic range where nonlinear effects are several 
times larger than at HERA. 

Electromagnetic interactions can also be studied with 
beams of protons or anti-protons, but there is then no 
Z^-enhancement in the photon flux in contrast to AA 
collisions. Several analysis are currently being done 
at Tevatron focusing on such processes. For instance, 
CDF Collaboration is analyzing the exclusive produc- 



tion of muon pairs, pp pp + n^p" , at lower invari- 
ant masses Q. The two main contributions to these 
events are, as with heavy-ion beams, 77 ^ ^+ -I- p~ 
and 7 + IP ^ J/^'(or4''), followed by the meson de- 
cay into a dilepton pair. The continuum 7 -I- 7 — > fJ-'^ fJ.~ 
has been calculated ^] under the same conditions as the 
vector mesons photoproduction and it was found a cross 
section, pp —> pp -\- fJ.~^ iJ-~ = 2.4 nb. In that calculation, 
it was used the cut Mjnv > 1-5 GeV for the dilepton 
invariant mass. 

In this work, we use the high energy color dipole ap- 
proach [5] to study the exclusive photoproduction of lep- 
ton pairs, 7iV -> 7*(-> £+£-) N (with N = p, A). The 
main reason is that such a process should be a back- 
ground for the main reaction 77 We use simple 
models for the elementary dipole-hadron scattering am- 
plitude that captures main features of the dependence on 
atomic number A, on energy and on momentum transfer 
t. This investigation is somewhat complementary to con- 
ventional partonic description of timelike Compton scat- 
tering (TCS) 0, 7P ^ 7* +p, which considers the rele- 
vant generalized parton distributions (GPDs). It should 
be noticed that the TCS process has so far only been 
studied at LO in the coUinear factorization framework 
in terms of the quark GPDs and sub-processes initiated 
by gluons have not been considered [6] . The color dipole 
approach provides a very good description of the data on 
7P inclusive production, 77 processes, diffractive deep in- 
elastic and vector meson production (For a recent review, 
see 7]). In particular, the electromagnetic deeply virtual 
Compton scattering (DVCS) cross section is nicely re- 
produced in several irnplementations of the dipole cross 
section at low x [sl, i, [l^, Ull, HI [11 • The TCS process 
at small t and large timelike virtuality of the outgoing 
photon shares many features of DVCS. 

This paper is organized as follows. In next section, 
we make a short summary on the color dipole approach 
applied to the exclusive protoproduction of dileptons at 
high energies. In Sec. 3, wc present the main results 
for such a process using proton and nucleus targets. For 
sake of simplicity in numerical calculations, we consider 
the space-like kinematics in this preliminary estimation 
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of relevant cross sections. Finally, in last section we sum- 
marize the results. 



II. TCS PROCESS WITHIN THE COLOR 
DIPOLE APPROACH 



Let us summarize the relevant formulas in the color 
dipole picture for the TCS process. In this formalism 
the scattering process 7p j*p is assumed to proceed 
in three stages: first the incoming real photon fluctuates 
into a quark-antiquark pair, then the qq pair scatters 
elastically on the proton, and finally the qq pair recom- 
bines to form a virtual photon (which subsequently de- 
cays into lepton pairs) . The amplitude for production of 
the exclusive final state such as a virtual photon in TCS, 
is given by 

f h.h-' ^0 
X Aqq{x,r,A)-^^-^{r,z,{)), (1) 

where f^f^{r, z,Q) denotes the amplitude for a photon 
to fluctuate into a quark-antiquark dipole with helicities 
h and h and flavour /. The quantity Aqq{x,r, IS) is the 
elementary amplitude for the scattering of a dipole of size 
r on the proton, A denotes the transverse momentum 
lost by the outgoing proton (with t = — A^), x is the 
Bjorken variable and is the photon virtuality. The 
elementary elastic amplitude Aqq can be related to the 
S'-matrix element 5'(x, r, h) for the scattering of a dipole 
of size r at impact parameter b [l^, [lH : 



Aqq{x,r, A) 



j A^b e-'^-^ 2[1 - S{x,r, b)] 



(2) 



As one has a real photon at the initial state, only the 
transversely polarized overlap function contributes to the 
cross section. Summed over the quark helicities, for a 
given quark flavour / it is given by (isj , 
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m'jKo{eir)Kois2r)} , 



(3) 



where we have defined the quantities 2 = zz Q\ 2 
and z — (1 — z). Accordingly, the photon virtualities are 
Q\ = Q (incoming real photon) and Q2 — —Q^ (outgo- 
ing virtual photon). For sake of simplicity, in numerical 
calculations we use space-like kinematics. We expect not 
a large deviation from the correct kinematics. However, 
the approximation we have considered to estimate the 
cross section should be taken with due care. We quote 
Ref. [ll] for a detailed derivation of wavefunctions for 
time-like photons. The elastic diffractive cross section is 
then given by (disregarding the real part of amplitude) , 



dt ~ 16^ 



^TP-'^'P(x,Q,A) 



(4) 



It should be noticed that some corrections to this ex- 
clusive process are needed. For TCS one should use the 
ofii'-diagonal gluon distribution, since the exchanged glu- 
ons carry different fractions x and x' of the proton's 
(light-cone) momentum. The skewed effect can be ac- 
counted for, in the limit that x' ^ a; ^ 1, by multiplying 
the elastic differential cross section by a factor i?^, given 
by [14] 
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In our numerical calculations we take into account sat- 
uration models which successfully describe exclusive pro- 
cesses at high energies (vector meson production, diffrac- 
tive DIS and DVCS). As a baseline model we consider 
the non-forward saturation model of Ref. [ll[ (hereafter 
MPS model). It has the great advantage of giving di- 
rectly the t dependence of elastic differential cross sec- 
tion without the necessity of considerations about the 
impact parameter details of the process. Furthermore, 
the overall normalization of the total cross section is com- 
pletely determined and no assumption about elastic slope 
is needed. Therefore, this simple saturation model cap- 
tures main features of the dependence on energy, vir- 
tual photon virtuality and momentum transfer t. An 
important result about the growth of the dipole ampli- 
tude towards the saturation regime is the geometric scal- 
ing regime [l^, [lB| ■ It first appeared in the context of the 
proton structure function, which involves the dipole scat- 
tering amplitude at zero momentum transfer. At small 
values of x, instead of being a function of both the vari- 
ables r and a;, the dipole scattering amplitude is actually 
a function of the single variable r^Q^^^^ix) up to inverse 
dipole sizes significantly larger than the saturation scale 
(5sat(a;)- More precisely, one can write 



Aqq(x, r, A = 0) = 2^i?2 N {r-'Ql^,(x)) , 



(5) 



implying the geometric scaling of the total cross-section 

at smafl x, i.e. alJ~'^(x,Q'^) = cTtV^^lTp^QV^D- 
This geometric scaling property can be extended to the 
case of non zero momentum transfer (l7| , provided r A ^ 
1. Therefore, it has been obtained that equation ^ can 
be generalized to 

Aqq{x,T,l\) ^2'KRle-''\'\^^ {tQ,^,{xM),^) . (6) 

with the asymptotic behaviors (5sa^t(x, A) ~ 
max((5Q, A^) exp[— Aln(a:)]. Specifically, the t de- 
pendence of the saturation scale is parametrised as 



' sat 



{x,\t\) = Qi(\^c\t\) ( - 



(7) 



in order to interpolate smoothly between the small and 
intermediate transfer regions. The form factor -F(A) = 
exp(— i3|t|) catches the transfer dependence of the pro- 
ton vertex, which is factorised from the projectile vertices 
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and does not spoil the geometric scaling properties. Fi- 
nally, the scaling function N is obtained from the forward 
saturation model [ll] , whose functional form is given by, 



N{x, r) 



1 _ Q-A\J{BrQ,^^^) 



2(7. + ^ In ^) 



rQsat < 2 
^Qsat > 2 

(8) 



where Y = ln(l/a;). The parameters for N are taken 
from the original CGC model [l8l |. 



III. PHENOMENOLOGICAL RESULTS FOR 
TCS AT HIGH ENERGIES 
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In this section the numerical results for the TCS cross 
section are presented using the color dipole picture and 
using the MPS saturation model as a representative 
phenomenological model. The sensitivity to a differ- 
ent choice for the dipole cross section is also analyzed. 
Here, we investigate the exclusive photoproduction of a 
heavy timelike photon which decays into a Icpton pair, 
7P — > £^£^p. Therefore, for the invariant mass 

distribution from the virtual 7* decay we have (with 

Q2 = M,V), 



da 



e-p) - 



. cr ^ 7» -(9) 



In Fig. [T] the differential cross section, (Pa/d\t\dA''P , 
is presented as a function the momentum transfer (with 
\t\ < 1) for representative values of the dilepton invari- 
ant mass. We consider the typical DESY-HERA en- 
ergy, Wjp — 82 GeV, where the DVCS cross section is 
currently measured. From a qualitative point of view, 
at small \t\ the differential cross section is consistent 
with an exponential behavior and a M^-dependent slope. 
The general \t\ dependence is determined by the ele- 
mentary dipole-nucleon scattering amplitude. For small 
size dipoles, one can qualitatively compute the scattering 
amphtude, which gives A oc (Q^^^/M^). For the non- 
forward saturation model, Qg^t oc (1 -I- c|i|) exp(— 
Thus, we checked that the cross section can be parame- 
terized as da/dtdM"^ = a{l + c\t\) exp(— 6|i|), with the 
constants a, c and b depending logarithmically on Af^, 
b = h{M'^). It should be noticed that the |t|-dependence 
is distinct for other implementations of the dipole cross 
section. However, for integrated cross sections, these de- 
viations probably do not play an important role. 

In the following we compute the integrated cross sec- 
tion, performing the integration over |i| < 1 GeV^. In 
Fig. [2] we focus on the invariant mass dependence for 
fixed values of energy. In order to test the energy range 
which can be covered in electromagnetic interactions in 
pp collisions, we start with the typical DESY-HERA 
value, VF = 82 GeV, and extrapolate it up to = 10 
TeV. Accordingly, the spectrum is dominated by small 



FIG. 1: (Color online) The differential cross section 
as a function of \t\ at fixed values of dilepton 
invariant mass. 
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FIG. 2; (Color online) The differential cross section, da/dM^ , 
as a function of dilepton invariant mass for fixed values of 
energy (integrated over |t| < 1 GeV^). 



invariant masses of order < 20 GeV'^, indepen- 

dent of energy. In the saturation models, the qualitative 
behavior can be obtained considering dominance of small 
size dipoles configuration. The forward amplitude reads 
as .4 cx (Qsat/-^^) times a logarithmic enhancement on 
NP . Therefore, qualitatively the differential cross sec- 
tion behaves like da/dM"^ cx {M'^)-^[l + \og{M'^)], with 
(5« 3. 

In Fig. [3] we compute the total cross section, cr(7p 
p), integrated over dilepton invariant mass Mg+g- > 
1.5 GeV. In this plot, we compare the MPS model 
(solid line) with other implementations of the elemen- 
tary dipole-nucleon scattering amplitude. For energies 
W >> 100 GeV, a power fit can be performed in the 
form a = A{W/Wo)°', with Wo = 1 GeV. Considering 
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the MPS model, one obtains the values A = 3 ph and 
a = 0.46. In order to study the sensitivity to the model 
dependence, we compare the MPS saturation model to 
two different distinct saturation models. The first one 
is the recent implementation of the impact parameter 
Color Glass Condensate model jTl] (hereafter WATT, 
dot-dashed line). The second one is the impact parameter 
saturation model [lOj Thereafter b-SAT, long dashed line). 
In these implementations, the elementary dipole-nucleon 
scattering amplitude is written in the impact parameter 
space as referred in Eq. In particular, in the b-SAT 
model the ^-matrix element is given by [10|, 



S{x, r, b) — exp 



2N, 



(10) 



Here, the scale is related to the dipole size r by 
^'^ = i/r'^ + ^Q. The gluon density, xg{x,y?), is evolved 
from a scale /ig up to using LO DGLAP evolution 
without quarks. The initial gluon density at the scale 
/^g is taken in the form xg{x,ijQ) = AgX^'^s (1 — x)^-^. 
The proton shape function T(b) is normalized so that 
Jd^b T{b) — 1 and one considers a Gaussian form for 

T{b), that is, Taib) = 2;^^^ where Bg = 4 GeV"^ 

At high energies, a power fit can be also performed for 
the WATT and b-SAT models. For WATT impact pa- 
rameter model, we obtain the parameters A = 4.3 pb and 
a = 0.37. In the b-SAT model we have the values A = 2 
pb and a = 0.57. We clearly verify a distinct energy de- 
pendence, which depends on the characteristics features 
of the phenomenological models. The main point is the 
value of the parameter A entering at the saturation scale. 
The WATT model gives the softer energy dependence, 
which comes from the small value of A = 0.12 in the sat- 
uration scale, Qgjjt(x,t = 0) oc a;'^. On the other hand, 
in the MPS model one has A — 0.22 (a factor two larger 
than in WATT model) and the QCD evolution makes 
the effective A value for b-SAT model to be large. The 
different overall normalization at energies W < 100 GeV 
for the MPS model is probably result of distinct behavior 
towards low energies. 

Let us discuss the connection of the present calcula- 
tion with the conventional partonic description of time- 
like Compton scattering. Presently, the photoproduc- 
tion of a heavy timelike photon which decays into a lep- 
ton pair, 7p £^£^ p, is computed to leading twist 
and at Born level 6]. In the Bjorken limit, i.e. large 
Q^, the QCD factorization theorem represents the am- 
plitude by the convolution of hard scattering coefficients, 
calculable in perturbation theory, and generalized par- 
ton distributions (GPDs), which describe the nonper- 
turbative physics of the process. GPDs are universal, 
process-independent, functions that contain information 
on parton distributions and correlations in hadrons and 
in matrix elements describing transitions between differ- 
ent hadrons. For instance, in Ref. [gJ, the quark hand- 
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FIG. 3: (Color online) The TCS integrated cross section 
{M(+i- > 1.5 GeV and \t\ < 1 GeV^) as a function of photon- 
nucleus centre of mass energy (see text). 



bag diagrams (leading order in as) and simple models of 
the relevant GPDs are used to estimate the cross section 
and the angular asymmetries for lepton pair production. 
We remark that such a calculation can not be compared 
with ours since an explicit analysis using 0{as) accu- 
racy, one-loop corrections to the quark handbag diagrams 
and other diagrams involving gluon distributions is cur- 
rently unknown. Finally, it should be noticed that the 
Bethe-Heitler (BH) mechanism contributes at the ampli- 
tude level to the physical process of photoproduction of 
heavy lepton pair. It is known that the BH contribu- 
tion (and its interference with TCS) is large in contrast 
to timelike Compton scattering and a careful numerical 
analysis is deserved. 

Finally, we present a study on the TCS process for 
scattering off nucleus. Namely, we focus on the process 
J A 7* {-^ £~^£~) A. In this case, we are considering the 
coherent contribution to nuclear TCS, where the nucleus 
remains intact. The implementation of nuclear effects in 
such a process is relatively simple within the color dipole 
formalism in the small- a; region. The usual procedure 
is to consider the Glauber-Gribov formalism for nuclear 
absorption. In phenomenological models in which geo- 
metric scaling is present (as in MPS and WATT) the 
extrapolation for a nucleus target is simplified. There- 
fore, here we will rely on the geometric scaling prop- 
erty [1^, [1^ of the saturation models: such a scaling 
means that the total 7*p cross section at large ener- 
gies is not a function of the two independent variables 
X and Q, but is rather a function of the single vari- 
able Tp = Q^/Q-^^^{x) as shown in Ref. [l^. That is, 
aj*p{x,Q'^) — crj*p{Tp). In Refs. [l^ it was shown that 
the geometric scaling observed in experimental data can 
be understood theoretically in the context of non-linear 
QCD evolution with fixed and running coupling. Re- 
cently, the high energy l^p, pA and AA collisions have 
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been related through geometric scaUng [l9|. Within the 
color dipole picture and making use of a rescaling of the 
impact parameter of the 7*/i cross section in terms of 
hadronic target radius Rh, the nuclear dependence of 
the 7* A cross section is absorbed in the A-dependence 
of the saturation scale via geometric scaling. The rela- 
tion reads as cr7o/(a;,Q^) = HAcr^f {Q sat, p Qsat,A), 
where ka = {Ra/ Rp)^ ■ The nuclear saturation scale was 
assumed to rise with the quotient of the transverse par- 
ton densities to the power Awl and Ra is the nuclear 
radius, Qg^t A ~ i^/ '^a)'^ Qsat,p- This assumption suc- 
cessfully describes small- a; data for ep and eA scattering 
using A = 1.26 and a same scaling curve for the proton 
and nucleus UM- 



Following the arguments above, we will replace 
Rp Ra in Eq. ^ and also Ql^^p{x,t = 0) ^ 
{ARl/R\)^ Ql^t.p{x, t = 0) in Eq. 0. In case of A = 1 
the previous replacement becomes the usual assumption 
for the nuclear saturation scale, QLt A — Qsat,p- For 
simphcity, we replace the form factor F{t) — exp(— 
in Eq. ^ by the corresponding nuclear form factor 

FA{t) = exp( — g^|t|). In a careful calculation, a real- 
istic nuclear form factor (for instance, hard sphere pro- 
file) is needed. The corresponding phenomenological re- 
sults are presented in Fig. [U where we have the in- 
tegrated cross section {Mi+i- > 1.5 GeV and \t\ < 1 
GeV^) per nucleon as a function of energy. We consider 
A = 208 (Lead), which is relevant for electromagnetic 
interactions at A A collisions at the LHC. The nuclear 
version of MPS model (using geometric scaling prop- 
erty) is represented by the solid line. At high energies 
{W > 100 GeV) the numerical result can be param- 
eterized as aupsilA A) = 6.1ph{W/Wof-^^. 
For the b-SAT model, we just replace the proton shape 
function T{b) in Eq (jTU]) by the corresponding nuclear 
profile TA{b) (Wood-Saxon). The result is presented by 
dot-dashed line and at high energies is parametrized as 
crbSAT(7A ^ e+£-A) = 5pb (W^/Wo)°-^^ Accordingly, 
we see a larger suppression in the MPS model than in 
b-SAT, which is directly related to the nuclear satura- 
tion scale at each model. Absorption is evident in MPS 
model, where the effective power on energy has dimin- 
ished in the nuclear case. 

Here, some comments are in order. The approximation 
we have considered to compute the photonuclear cross 
section should be taken with due care. In the small-x 
limit, within the GPDs formalism the ratio of the real 
parts of the nuclear to nucleon amplitudes exclusive pro- 
cesses (like DVCS) has a very unexpected behavior as 
a function of x, which is distinct from the correspond- 
ing ratio of imaginary parts. This is the reason we have 
disregarded the real part contribution in our numerical 
calculations. Another important aspect is the interplay 
between the coherent and incoherent contributions to the 
nuclear TCS. When the recoiled nucleus is not detected, 
the measurements of TCS observables with nuclear tar- 
gets involves both contributions and their interference. 



In the coherent scattering the nuclear target stays intact 
and it dominates at small-t. On the other hand, in the 
incoherent process the initial nucleus A transforms into 
the system oi A—1 spectator nucleon and one interacting 
nucleon (it dominates at large-i). We have considered in 
present study the situation when the recoiled nucleus is 
detected and the 7 A ^ ^*A cross section is purely co- 
herent. For a scenario with excited or break up nucleus, 
the cut \t\ ^1 GeV^ could prevent us of considering 
the incoherent contribution. The coherent contribution 
dominates the nuclear cross section at small-t. 

The present calculations are input for the exclusive 
photoproduction of dileptons in electromagnetic inter- 
actions in pp and nucleus- nucleus collisions. Such pro- 
cesses are characterized by the photon - hadron inter- 
action, with the photon stemming from the electromag- 
netic field of one of the two colliding hadrons (For recent 
reviews see Ref. [l[). The total cross section for the 
hh ^ h (S) £~^£~ (g) h process is given by 

ihh - £^£- hh) ^ 2^-^ (7/. -> £n- h) 



did dM^^^_ 



where uj is the photon energy and — is the equivalent 
flux of photons from a charged hadron. Moreover, 7^^ is 

2 wV'S'nn and 



the Lorentz boost of a single beam, W? 



7/1 

vSnn is the c.m.s energy of the hadron-hadron system 
n|. This process is characterized by small momentum 
transfer and energy loss, which implies that the outgo- 
ing hadrons should be detected in the forward regions 
of detectors. A rough estimation for pp collisions gives 
dozens of pb at Tevatron and a hundred of pb at the 
LHC (integrated over invariant mass Mi^v > 1-5 GeV). 
In the AA the photon flux is enhanced by a factor (x Z^, 
then we roughly expect the cross sections reach dozens of 
nb at PbPb collisions at the LHC. The exclusive photo- 
production of dileptons from TCS has similar final state 
configuration than the QED process 77 — > How- 
ever, detailed estimations deserve further studies as the 
correct time- like kinematics and realistic nuclear profiles, 
which we postpone for a future publication. 



IV. SUMMARY 

Using the color dipole formalism, we studied the time- 
like Compton scattering. Such an approach is robust in 
describing a wide class of exclusive processes measured 
at DESY-HERA and at the experiment CLAS (Jeferson 
Lab.), like meson production, diffractive DIS and DVCS. 
Our investigation is complementary to conventional par- 
tonic description of TCS, which considers quark hand- 
bag diagrams (leading order in as) and simple models of 
the relevant GPDs. In particular, the results could be 
compared to pQCD diagrams involving gluon distribu- 
tions which are currently unknown. Using current phe- 
nomenology for the elementary dipole-hadron scattering, 
we estimate the order of magnitude of the exclusive pho- 
toproduction of lepton pairs. In the present analysis we 
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M,„„>1.5 GeV 
A = Pb 



MPS (nuclear) 
b-SAT (nuclear) 



that the cross section per nucleon reaches dozens of pi- 
cobarns. At DESY-HERA (with W ^ %2 GeV) the pure 
TCS contribution gives a cross section of 20 pb with an 
uncertainty of a few percents. It should be noticed that 
we consider space-hke kinematics for numerical calcula- 
tions and that the correct time- like kinematics will intro- 
duce corrections to our preliminary estimations. These 
calculations are input in electromagnetic interactions in 
pp and AA collisions to measured at the LHC. We found 
that the exclusive photoproduction of lepton pairs in such 
reactions should be sizable. 
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FIG. 4: (Color online) The photonuclear cross section per 
nucleon as a function of energy for a lead nucleus (see text). 

also investigate the photonuclear cross section, focusing 
on the nuclear coherent scattering. The main results is 
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